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. ABSTRACT 

' We investigate the science goals achievable with the upcoming generation of ground- 

^ I based Cosmic Microwave Background polarization experiments, focusing on one partic- 

^ ular experiment, QUaD, a proposed bolometric polarimeter operating from the South 

Pole. We calculate the optimal sky coverage for this experiment including the effects 
of foregrounds and gravitational lensing. We find that an E-mode measurement will 



o 



X 



be sample-limited, while a B-mode measurement will be detector-noise-limited. We 
^ I conclude that a 300 deg^ survey is an optimal compromise for a two-year experiment 

to measure both E and B-niodes, and that a ground-based polarization experiment 
can make an important contribution to B-mode surveys. QUaD can make a high sig- 
nificance measurement of the acoustic peaks in the E-mode spectrum, over a multipole 
range of 25 < £ < 2500, and will be able to detect the gravitational lensing signal in 
' the B-mode spectrum. Such an experiment could also directly detect the gravitational 

wave component of the B-mode spectrum if the amplitude of the signal is close to 
current upper limits. We also investigate how QUaD can improve constraints on the 
cosmological parameters. We estimate that combining two years of QUaD data with 
the four-year WMAP data can improve constraints on ftbh^, flmh^, h, r and Ug by 
a factor of two. If the foreground contamination can be reduced, the measurement of 
r can be improved by up to a factor of six over that obtainable from WMAP alone. 
These improved accuracies will place strong constraints on the potential of the inflaton 
field. 

Key words: cosmic microwave background - cosmological parameters - 
methods:observational - polarization - techniques ipolarmetric 
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1 INTRODUCTION 

The Cosmic Microwave Background (CMB) has proven to 
be a powerful cosmological probe. Successive generations of 
experiments have provided a stringent test for the standard 
Big Bang paradigm and increasingly sensitive measurements 
of the temperature power anisotropics have led to tight con- 
straints on many of the fundamental cosmological parame- 
ters. However, as well as fluctuations in the CMB tempera- 
ture field, there are also anisotropies in the linear polariza- 
tion of the CMB. These polarization fluctuations have re- 
cently been detected by the DASI experiment (Kovac et al. 
2002) and the correlation between the temperature and po- 
larization has been measured by the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite (Kogut et al. 2003). 
However, to make full use of the CMB, higher sensitivity 
high resolution polarized measurements arc needed. This is 
the challenge facing the next generation of CMB experi- 
ments. 

While it is desirable to observe the CMB tempera- 
ture field from space, to remove atmospheric noise, this 
is not as important for polarization experiments since the 
atmospheric emission is not expected to be linearly polar- 
ized (Keating et al. 1998). Therefore, by integrating deeply 
on relatively small patches of sky (Jaffe, Kamionkowski & 
Wang, 2000) it is possible to make a measurement of the 
polarization anisotropies with a comparable signal-to-noise 
ratio to a satellite experiment on all but the largest angular 
scales. 

The survey design for a ground-based experiment will 
depend upon the specific science goals of the experiment. In 
this paper we investigate observing strategies and sky cov- 
erage for the forthcoming generation of ground-based CMB 
polarization experiments, taking into account foreground is- 
sues. We will also show how ground-based polarization mea- 
surements can help to tighten constraints on the cosmolog- 
ical model. 

To be concrete, we focus on one particular experiment, 
QUaD (QUEST: Q and U Extra-galactic Sub-millimetre 

Telescope, and DASI: Degree Angular Scale Interferometer). 

This is a proposal to install QUEST§ , a high-resolution bolo- 
metric array polarimeter, on the azimuth-elevation mount of 
the DASI^ instrument. The experiment plans to begin ob- 
serving from the South Pole in 2005 (Church et al. 2003). 

The remainder of the paper is set out as follows. In Sec- 
tion 2 we briefly review the physics of the CMB polarization. 
In Section 3 we present the formalism used in the investi- 
gation and in Section 4 we show how we have included the 
effects of foregrounds. Our cosmological model and defini- 
tions are presented in Section 5. In Section 6 we present our 
results for the survey design, and in Section 7 simulate po- 
larization maps. The expected accuracies and multipole cov- 
erage of the power spectra are presented in Section 8, and 
in Section 9 we present the expected parameters constraints 
for the QUaD experiment. Our findings are summarized in 
Section 10. We also include an appendix in which we discuss 
the sensitivity definitions used in our calculations. We begin 
with a brief review of CMB polarization. 



2 REVIEW OF CMB POLARIZATION 

Detailed reviews of the CMB polarization are given by Zal- 
darriaga (2003) and Hu & White (1997). In this Section we 
give a brief overview of how the polarization field is gener- 
ated and how it is parameterized. 



2.1 Parameterization of the polarization field 

Typically, a linearly polarized source is quantified by the Q 
and U Stokes parameters, which give the differences in in- 
tensity between orthogonal polarization states. These quan- 
tities are convenient to measure experimentally, but are not 
invariant under a rotation of the co-ordinate system and so 
are difficult to compare to theoretical models. It is useful 
to define the quantities Q(9,(f)) ± iU{6,(f)) which under ro- 
tation acts as a spin 2 quantity. By operating on Q ± ill 
using spin raising and lowering operaters it is possible to 
obtain two rotationally invariant spin fields, E{6, </>) and 
B{9,cf)) (Zaldarriaga & Seljak, 1997). This represents the 
decompostion of the polarization field into different parity 
states, the E field is unchanged by a parity transformation, 
but the B field changes sign. These E and B fields can then 
be compared directly to theoretical predictions H. 

The temperature field is usually expanded in terms of 
scalar spherical harmonics, Yim{0,4')'- 

em 

where AT is the deviation of the temperature field from its 
average value To. Similarly, the quantities Q ± ill can be 
expanded in terms of spin 2 spherical harmonics, :f2Yem- 

Q{e, <j>) ± iU{e, 0) = ^ {Eem T iBem)^2 ^(mie, 4>)- (2) 
em 

The two point statistics of the CMB can be completely 
described in terms of the covariances of the multipole mo- 
ments, Tim, Elm and Bim'- 

^ImTfrn') = Cj^ du'^mm' {E(^Eii^i) = Cf^ Su'^mm' 
{Blm^i'm') = C!( Su'Smm' (^/m-^^'m') = ^u'^mm' 
{Tem^e'Tn') = Ce Su>5mm' {Eljn^t'm') = Cg^ fe'<^mm'- 

As the B field has opposite parity to the T and E fields the 
TB and EB correlations are zero if we can assume that par- 
ity is conserved. If the CMB is a Gaussian random field, as 
predicted if the metric fluctuations are generated by zero- 
point fluctuations during inflation, the statistical properties 
of the CMB temperature and polarization fields are com- 
pletely defined by the four power spectra, Cj^ ,Cf^ ,Cf^ 
and Cj^ . However, as we shall discuss in Section 2.3, grav- 
itational lensing by large-scale structure along the line of 
sight will distort the pattern of fluctuations, and will induce 
non-Gaussianity. 



§ http://www.astro.cf.ac.uk/groups/instrumentation/projects/ 
^ http://astro.uchicago.edu/dasi/ 



II An equivalent formalism is given by Kamionkowski, Kosowsky 
& Stebbins (1997) in which the polarization field is expanded in 
terms of tensor spherical harmonics instead of spin 2 harmonics. 
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2.2 Polarization signal generated during 
recombination 

The CMB polarization signal primarily arises from the 
Thomson scattering of the CMB photons during recombi- 
nation. Polarization can only be generated if the radiation 
field contains a local quadrupolo. Density perturbations will 
produce a velocity gradient in the primordial plasma so that 
photons approaching an electron from different directions 
will be Doppler shifted by different amounts. This produces 
local quadrupoles in the radiation field. Before recombina- 
tion, the high electron density means that the mean free 
path of the photons is too small to produce a quadrupole; 
however, after the recombination the electron density is too 
low for significant Thomson scattering to occur. The polar- 
ization can only be produced during a short period around 
recombination, so the amplitude of the polarization is very 
low. 

The mechanism by which these scalar perturbations are 

produced in the polarization field is therefore subtly differ- 
ent to the way in which the temperature perturbations are 
produced. A measurement of the polarization power spectra 
will not only provide a consistency check of the cosmologi- 
cal model, but will also yield new information on processes 
occurring in the early universe. Much of this information is 
contained in the TE and EE acoustic peaks at high I which 
can be measured with high signal to noise with a ground- 
based experiment. 

The inflationary model also predicts a stochastic back- 
ground of gravitational waves (GW) which will also result 
in a quadrupolo. The decomposition of the polarization field 
into the E and B modes can be used to separate the GW 
(tensor) contribution from the density perturbation (scalar) 
contribution. The E-modes can be produced by both scalar 
and tensor perturbations, but the B-modcs produced at last 
scattering can only be generated by tensor perturbations. 
This means that a measurement of the B-mode spectrum 
would give new information about inflationary parameters. 
In particular, the amplitude of the tensor spectrum is di- 
rectly related to the energy scale of inflation. These param- 
eters can not be well constrained from the TT and EE spec- 
trum as it difficult to separate the tensor and scalar con- 
tributions to these measurements. The GW B-mode signal 
peaks around scales of about £ = 100 and so in principle is 
detectable from the ground. 

2.3 Polarization signal generated after 
recombination 

The polarization spectra generated at recombination will be 
altered mainly by two processes before they can be detected: 
re- ionization and weak gravitational lensing (GL). The effect 
of re-ionization is to increase the polarization signal on large 
scales (£ < 20). Ground-based experiments are unlikely to be 
able to measure the polarization on such large angular scales 
and so will not be sensitive to the effects of re-ionization. 
However, weak lensing affects the signal on small angular 
scales. CMB photons are deflected by the gravitational po- 
tential of large scale structure. For the TT and EE spectra, 
this effect results in a smearing of the acoustic peaks on small 
angular scales, although the change to the spectra is very 
small, as shown in Fig. 1. However, lensing will also convert 
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Figure 1. CMB temperature and polarization power spcetra. 
The dashed lines are for a model with no re-ionization while the 
dotted lines are for a model with no gravitational lensing. The 
solid lines include the effects of both gravitational lensing and 
re-ionization. Model parameters are given in Section 5. 



E-mode polarization into B-modes. This means that there 
will be a scalar contribution to the B-mode spectrum due 
to lensing. Therefore, tlic B-mode spectrum will be contam- 
inated by a GL contribution, the spectrum of which must 
be measured precisely so that it can be removed (Knox & 
Song 2002, Kesden, Cooray & Kamionkowski 2002). As the 
lensing signal peaks at small angular scales, ground-based 
experiments are well-suited to this task. 

The lensing signal itself also contains useful informa- 
tion about large-scale structure. This can be used to con- 
strain other cosmological parameters such as the neutrino 
mass (Kaplinghat, Knox & Song 2003), since this will add 
to the mass-energy of the universe, altering its expansion 
history, and suppressing small scale power in the matter 
power spectrum due to free streaming. The lensing signal 
will also make the CMB sensitive to the equation of state of 
the universe, parameterized by m; = pj p, as again this will 
affect the expansion history. 

Fig. 1 shows the temperature and polarization power 
spectra, generated by the Boltznianii and Einstein solver 
CMBFAST (version 4.2; Seljak & Zaldarriaga (1996)**, 
decomposed into temperature-temperature (TT) power, 
temperature-E-mode (TE) cross-power, E-mode-E-mode 
(EE) power, and B-mode- B-mode (BB) power. We plot spec- 
tra without gravitational lensing (dotted lines) and without 
re-ionization (dashed lines) and with both included (solid 
lines). The main aim of this paper is to determine the best 
survey design for the measurement of the polarization spec- 
tra. In the following section we present our formalism for 
this procedure, based on the Fisher Information matrix. 



http: / /www.cmbfast.org/ 
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3 FORMALISM 

3.1 Fisher Information matrix 

For a model dependent on a set of parameters, a, the proba- 
bility of a particular parameter set, given a set of experimen- 
tal data points, d, is expressed by the likelihood function, 
L(a|d), the probability of the parameters given the data. 
By exploring the parameter space to maximize L wo may 
determine the parameter values within certain error limits. 
The minimum possible variance with which a parameter can 
be measured can be estimated from the Fisher information 
matrix (Tegmark, Taylor & Heavens 1997), defined as: 



^ daidai 



(4) 



where C = — In L and the derivatives arc evaluated at the 
maximum likelihood values of the parameters. The inverse 
of the Fisher matrix gives the parameter covariance matrix, 
Cij, for the theoretical parameters: 



Cij = (AaiAaj) ^T, 



-1 



(5) 



where Aa^ is the deviation of the parameter from its max- 
imum likelihood value. The diagonal of the inverse Fisher 

matrix yields the marginalized \-a error on the parameters. 
Taking the inverse of the diagonal of the Fisher matrix, 

(AaO' = IjTu, (6) 

yields the conditional error on the parameters. In general 

> lITru (7) 

where the equality holds only for uncorrelated parameters. 
The Fisher matrix then provides a theoretical upper bound 
on the accuracy of a measurement of a given parameter for 
a given experiment. 



3.2 Application of Fisher matrix to CMB 
experiments 

For a CMB experiment, the data are the measurements of 
the four CMB power spectra and the parameters are the 

cosmological parameters. For the measurement of a single 
power spectrum, Ct, the Fisher matrix is given by: 



E 



1 



dCe dCe 



(AC£)2 dai daj ' 



(8) 



where 



i2£ + l)UyA£'- 

is the error in the measurement of the power spectrum in a 
band centred on multipolo £, and A'^ is a noise term. The 
survey area is given by fsky The summation is over pass- 
bands of width A£. 

For a measurement of all four power spectra this gen- 
eralizes to (Zaldarriaga & Seljak 1997): 



\XY 



del 

da-j 



(9) 



where X and Y are either TT, EE, TE or BB and Hxy = 
Cov{Cf CY) is the power spectra covariance matrix: 
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^TT,EE 


■^TT,TE 
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,^EE,EE 
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^TT,TE 


■;^EE,TE 
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~;BB,BB 



(10) 



The terms in the power spectra covariance matrix are given 
by: 



{2£ + l)UyA£ 



+{cr +ivr )(cr +^1 



(11) 



where {x,y) = {T,E,B). The noise covariance is given by 



N^'^ . In the case of no foregrounds this is given by: 



xy, (12) 

where, = ^pixi^pix)^ , for an experiment with solid an- 
gle per pixel, ^pix, and the noise per pixel, apix- The pixel 
noise depends on survey design and instrument parameters. 
For an experiment covering an area for an integration 
time tabs, with iVpge detectors, a solid angle per pixel iipix 
and a sensitivityl^^^, NET, the pixel noise is: 

ro6sJVpSB"pia; 

In Section 4 we discuss how the noise terms may be extended 
to include foregrounds. We assume that the pixel size used 
in the map wil be the same as the beam size of the telescope. 
The spherical harmonic transform of the beam is given by 
Be- Here we assume that the beam is a Gaussian, 



Bt = cyiY>{-({(+l)(jl/2) 



(14) 



with (JB = Ob/ V81n2 where 9b is the full width half max- 
imum beam size. The pixel size can then be approximated 

by Vlpix = Og. 

The minimum resolution of the power spectra, A£, de- 
pends on the area of sky covered, Al = n/Q. This will there- 
fore also give the minimum i at which the power spectra 
can be measured as discussed further in Section 6.2. If a 
resolution smaller than this is used, the different £ modes 
will become correlated and equation (3) will no longer ap- 
ply (Hobson & Magucijo 1996). We calculate the maximum 
£ value from the FWHM beam size, £max = ir/O. In reality, 
multipoles higher than this could be measured if the beam 
profiles can be accurately determined. 

The Fisher matrix also provides a simple way to cal- 
culate the results obtainable by combining a number of ob- 
servations from different CMB experiments. In the simplest 
case, in which Nexp experiments observe different patches 
of sky, the combined Fisher matrix, JT'-^, is the sum of the 
individual Fisher matrices, J^" (Hu 2001): 



(15) 



If any of the patches of sky overlap, each overlapping region 
is considered as a separate patch. In these patches the com- 
bined noise covariance, N^, of the overlapping experiments 

tt The definition of sensitivity for a polcirization experiment is 
discussed in Appendix (A). 
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should be used to calculate the terms in the power spectra 
covariaiice matrix (equation (11)). This is discussed further 
in the next section where we consider how to optimally com- 
bine multi-frequency data. 

This completes the formal machinery wc will require 
for our analysis. Note that we have ignored the effects of 
windowing and mode-mixing due to limited sky coverage 
(e.g. Bunn 2002), and non-Gaussianity and mode-coupling 
induced by gravitational lensing (e.g. Guzik, Seljak & Zal- 
darriaga 1999). The former effects modes by convolving 
them with the survey window function and mixing E and 
B modes. This will mainly effect the B-modes, where the 
signal-to-noisc is poor, and will slightly increase our uncer- 
tainties. Non-Gaussianity induced by gravitational lensing 
will also correlate modes and will give rise to higher-order 
correlations, which will also lead to a slight increase in our 
uncertainties. 

So far we have also ignored the effects of foreground 
contamination, and it is to this we now turn. 



4 FOREGROUNDS 

4.1 Including foregrounds into the formalism 

The signal measured from the sky will contain not only a 
component from the CMB, but also a contribution from 
astrophysical foregrounds. The CMB signal is independent 
of the wavelength of the observation, but the signal from 
most foregrounds is expected to be frequency dependent. 
By observing in a number of different frequency channels 
it is therefore possible to reduce the total foreground con- 
tamination by optimally combining the signal from different 
frequency channels. It may also be possible to use the multi- 
ple frequency information to remove some of the foreground 
contamination from the signal (e.g. Maino et al. 2002, Hob- 
son et al. 1998). 

The effect of observing over multiple channels needs to 
be taken into account in the Fisher matrix formalism de- 
scribed in the previous Section. If we ignore foregrounds and 
consider only detector noise, we can simply replace the noise 
terms in equation (11) by an inverse variance weighting of 
the noise in each channel, Nt c- 



(16) 



By choosing this weighting scheme at each multipole we 
combine the signals by giving the most weight to the chan- 
nels with the smallest detector noise. 

We include the effect of foregrounds by treating the fore- 
grounds as an extra source of noise with power spectra A?/® 
for each different power spectra in each frequency channel. 
This gives us the maximum possible foreground contamina- 
tion i.e. the contamination assuming that no foreground re- 
moval will be attempted. However, unlike the detector noise, 
the foregrounds will be correlated between power spectra 
and between frequency channels. To include these correla- 
tions we follow the technique developed in Tegmark et al. 
(2000, hereafter TOO). We define a 3F x 3F noise matrix. 
Mi, for each multipole, where F is the number of frequency 
channels in the experiment: 








XX' 



(17) 



is an F X F 



where each component of this matrix, N 
matrix giving the variances and covariances of the noise in 
the F channels. Each element in Mt is the sum of the con- 
tribution from each of the possible foregrounds, and 

the detector noise, N 



N 



XX' 
e(det) 



E 



XX' , 

e(det)- 

XX' 



N 



(18) 



where the sum over A: is a sum over each of possible fore- 
grounds which could contribute to the signal. We define the 
3F X 3 scan matrix. A, where: 



(19) 



and e is a column vector of height F consisting entirely of 
ones. If F = 2 as would be the case for QUaD (see Section 
6.2) then: 




A = 



( 



\ 



(20) 



The weighted noise for each polarization is then obtained by 
calculating the 3x3 covariance matrix, E^, where: 



Nr 




















Nr 



(21) 



The terms, are now the noise terms used in equation 

(11) to calculate the power spectra covariance matrix. If the 
noise is not correlated between T and E and not correlated 
between channels (as is the case if we include only detector 
noise) then Ne becomes diagonal and the procedure is iden- 
tical to the minimum variance weighting of equation (16). 

In the last section we discussed how to combine a num- 
ber of experiments by adding the Fisher matrices of inde- 
pendent patches of sky. For patches in which a number of 
experiments overlap, the required noise term, Ne, can be 
calculated by considering a single experiment with chan- 
nels at each of the different frequencies used by this set of 
experiments. For this patch F will then become the total 
number of frequency channels in the combined survey. If 
any of the instruments used have channels which cannot 
measure either temperature or polarization, then rows and 
columns corresponding to these channels should be removed 
from the full noise matrix, N and from the scan matrix, 
A, in the relevant places. For example, QUaD would not be 
able to measure temperature information (see Section 6.2). 
If we combine the two QUaD channels with another exper- 
iment measuring both temperature and polarization, these 
two channels should be removed from the first row and first 
column of the matrix A/i in equation (17) and the size of 
the vector e in the first column of the matrix A in equation 
(19) should be reduced. 
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Figure 2. Models used for vibrating dust, synchroton emission 
and residual point sources (black) compared to the EE and BB 
power spectra (light grey, solid). The different lines show fore- 
ground models at 100 GHz (dash-dot) and 150 GHz (dash). The 
total foreground power spectra are also shown on each plot (dark 
grey). The GW component of the BB-spectra is shown for r = 0.1 
(upper dotted line) and r = 0.01 (lower dotted line). The other 
foregrounds are either unpolarized or can be neglected 



4.2 Foreground models 

We closely follow TOO in constructing the foreground power 
spectra required in the previous section, and use the software 
provided on the associated website QUaD proposes to ob- 
serve at frequencies of 100 and 150 GHz. At these frequen- 
cies the relevant foregrounds are diffuse free- free emission, 
IR and radio point sources, synchrotron radiation and vi- 
brating dust, rotating dust and thermal Sunyaev-Zeldovich 
(SZ) radiation. Each foreground is modelled using a spatial 
power spectrum, Ci (k) = {pA)'^ , where f3 gives the scale 
dependence of the foreground fluctuations, p is the fraction 
polarized and A is the overall amplitude. A frequency de- 
pendence is also defined and normalized to unity at a refer- 
ence frequency, j/* . For point sources it is assumed that very 
bright sources (5 a outlyers) will be removed from the CMB 
maps, but that there will still be a residual point source con- 
tamination after this subtraction. In TOO, sets of estimates 
for these parameters are given. We begin by using their 
"middle-of-the-road" foreground model. In this model, the 
only polarized foregrounds are sychrotron, dust and point 
sources. We then slightly modify this model to take into ac- 
count recent observations (Kovac et al. 2002; Bennett et al. 
2003, hereafter BOS). These modifications lower the ampli- 
tude of the vibrating dust component and slightly increase 
the amplitude of the synchrotron emission. The amplitudes 
then roughly match those given in Fig. 10 of B03. Also fol- 
lowing B03 we have lowered the amplitude of the radio point 
sources and have neglected rotating dust emission. The val- 
ues of those foreground parameters which are different from 
TOO are given in Table 1. 

The power spectra of the relevant foreground models 



http:/ /www. hep. upenn.edu/^max/foregrounds. html 




Figure 3. Equal-area zenithal projection showing foreground lev- 
els (dust and synchrotron) at 150 GHz in regions which would be 
accessible to QUaD. The Southern Celestial Pole in located in 
the centre of the plot, and declination -45 around the perime- 
ter. To the bottom is Right Ascension 0, increasing in RA in the 
anti-clockwise direction. The possible QUaD observing region is 
shown by solid lined box. 



Table 1. Parameters used in foreground models which differ from 
those used in TOO. All other parameters used are as per the 
"middle-of-the-road" model in TOO. 



Foreground 


Radio point 


Synchrotron 


Vibrating Dust 




sources 






A/fJ.K 


0.66 


95 


7.5 


J/, /GHz 




20 


90 



are shown in Fig. 2 for the two QUaD frequency bands. The 
sychrotron radiation dominates the foregrounds at 100 GHz, 
whereas at 150 GHz both vibrating dust and sychroton ra- 
diation are important. The points sources only contribute 
at very high multipoles. For most of the multipole range of 
interest the EE spectrum dominates over the foregrounds. 
However, for the smaller BB signal the total foreground con- 
tamination is larger than the signal of interest. 

The analysis described in Section 4.1 gives the resid- 
ual foreground contamination given that foreground power 
spectra are well known or can be measured from the ex- 
perimental data. For QUaD we assume that this is rea- 
sonable given that other experiments, e.g. WMAP (B03), 
ARCHEOPS (Benoit et al. 2003) and the recent Boomerang 
fiight (Montroy et al. 2003), will soon provide polarized maps 
at CMB frequencies. Recent advances in foreground removal 
techniques (Baccigalupi 2003) indicate that it may be possi- 
ble to remove some of the foreground noise from the signal. 
The residual foreground contamination used here therefore 
gives an upper limit on that which can be expected in the 
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final cleaned maps, given that our foreground models are 
accurate. 

For a ground-based experiment it will also be possible 
to select preferentially regions of sky to observe in which 
the foreground fluctuations are small, and so the foreground 
noise can be reduced further. Fig. 3 shows the region of sky 
which would be accessible to QUaD from the South Pole and 
the estimated levels of foreground contamination across this 
area from dust (Finkbeiner, Davis & Schlogcl 1999) and syn- 
chroton (Giardino ot al. 2002) using the modified foreground 
models. A possible observing patch for QUaD is also shown 
in which the mean foreground amplitude is low. However, a 
more detailed analysis will be performed to choose the final 
observing patch with the lowest possible foreground vari- 
ance across the mutipole range of interest. For these reasons 
we perform the relevant calculations once using the full-sky 
foreground models descibed here, and again assuming that 
the foregrounds are negigible. 

We conclude that while our current understand of po- 
larized foregrounds is evolving, the expected level of fore- 
ground contamination in the EE power spectrum should not 
be significant. The GL component of the BB power spec- 
trum should also be detectable if patches of sky with low 
foreground variance can be targeted or if foreground removal 
techniques can be successfully implemented. This would also 
mean that the GW B-mode component should be measure- 
able if the tensor-to-scalar ratio is large. 



5 COSMOLOGICAL MODEL 

In order to calculate the terms in the power spectrum covari- 
ance matrix and the power spectrum derivative we require a 
model from which to calculate the CMB power spectra. This 
model is defined by two sets of parameters: the inflation- 
ary parameters which parameterize the initial perturbations 
causing the fluctuations in the CMB, and the cosmological 
parameters which determine how these initial perturbations 
are propagated into the observed CMB power spectra. Given 
a set of parameters, the CMB power spectra can then be 
calculated using a Boltzmann and Einstein solver. For this 
work we have used a slightly modifled version of CMBFAST 
vA.2. 

The initial scalar perturbations are parameterized by: 



A^(fe)=A^(fco)(^)' 



(22) 



where A^(A;) is the power spectra of TZ, the curvature per- 
turbation in the comoving gauge, and Us is the slope of the 
scalar power spectrum. The tensor perturbations are given 
by: 



A|(fe) = A|.(fco) 



(23) 



where Ay(A;) is the power spectra of gravitational waves 
from inflation and nt is the slope of the gravitational wave 
power spectrum. The amplitude terms are evaluated at the 
pivot wave number, fco = 0.05Mpc~^. To parameterize the 
initial perturbations we use three inflationary parameters: 
A, a constant of order unity which is proportional to the 
amplitude of the initial scalar perturbations, Us , the tilt of 
the power spectra of the initial scalar perturbations and r, 



the ratio of tensor to scalar perturbations. These are the 
parameters used in the analysis of the WMAP data (Spergel 
et al. 2003). We do not consider here the running of the 
spectral index, n'^. The exact relationship between A and 
A^(fc()) is derived in Verde ct al (2003; equation (32)). The 
tensor-to-scalar ratio is defined as: 



(24) 



Note that a number of different definitions arc used in the 
literature. The most common alternatives are to define r in 
terms of the Newtonian potential: 



A|(fco) 
A^(feo)' 



(25) 



so that 
quadrupoles: 



TQ = 



(5/3)^r, or in terms of the CMB radiation 



(26) 



The relation between r and rq depends on the cosmological 
parameters used in the model (Turner & White 1996). 

The cosmological parameters we shall consider are 
{Q.th^ , flmh^, h, r}, where h is the Hubble constant in units 
of lOOkms^^Mpc^^, Q.h is the energy density of baryons, 
Q.m is the total matter density and r is the optical depth to 
the last scattering surface. Again, these are the parameters 
chosen for the WMAP data analysis (Verde et al. 2003). The 
full set of parameters, and is then: 

{Hb/i^, fim/i^, h, T, Us, r, A} = 

{0.0224, 0.135, 0.71, 0.17, 0.93, 0.01, 0.83}. 

The values of these parameters are taken from the best-fit 
WMAP model (Spergel et al. 2003). Although we do not 
include n'^ in our analysis we note that these best-fit param- 
eters from a model which inculdes a non-zero n'g. For this 
parameter set the relation between the r and vq is approxi- 
mately rg Ri 2.8r and A|j.(feo) = 2.45 x lO"''. The values of 
these parameters are taken from the best fit WMAP model 
(Spergel et al. 2003) except for r, which cannot be well con- 
strained by this data set. The current upper limit on r is 
about 0.36 (Leach & Liddle 2003)§§. The lowest possible r 
which can be detected is of the order of 10~* (Knox & Song 
2002, Kesden, Cooray & Kamionkowski 2002) due to noise 
left over from the removal of the gravitational lensing signal 
from the B-mode spectrum. To refiect this range of possible 
values we perform the calculations, which have strong de- 
pendence on r, at two difFerent values, r = 0.01 and r = 0.1. 



6 SURVEY DESIGN 
6.1 Method 

The optimization of the survey area for a ground based mea- 
surement of the CMB polarization has been addressed previ- 
ously (Jaffe et al. 2000) in the context of making a detection. 



SS Note that our value differs from the value given in this refer- 
ence as we use a different value for the pivot wavenumber, fco 
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We extend this work by considering the criteria for a mea- 
surement of the polarization spectra including the effects of 
gravitational lensing and foregrounds. 

The main aim of a polarization experiment is to make 
measurements of the throe polarization power spectra, Cf^ , 
Cf^ and Cf^, with the highest possible precision. The error 
in the measurement of the power spectra is determined by 
two conflicting factors. For a fixed total observing time, the 
integration time per unit area (or pixel) is inversely propor- 
tional to the total area; a smaller map will therefore result 
in a lower pixel noise. However, for a smaller map there are 
fewer independent modes from which to measure each mul- 
tipole (i.e. the averaging in equation (3) will be made over 
fewer values of m) and so the sample variance will increase. 

To quantify these effects, we choose a single parameter 
for which to evaluate the Fisher matrix, , the amplitude 
of each power spectrum. For a single parameter the variance 
in the measurement of this parameter, (AA^)'^, is then given 
by ^/ F^xj^x ■ From equation (8) the error in is: 



Table 2. Expected QUaD instrument parameters 



1 



(ACf )2 {A^Y 



(27) 



where (AC^)^ for each power spectrum are given by the 
diagonal elements of the power spectrum covariance matrix 
in equation (10). We then define a figure of merit parameter 
as the signal to noise ratio in the measurement of each power 
spectrum, SNR, which is given by: 



E 



ACf 



(28) 



To find the optimal area for a measurement of each power 
spectrum with a spccfic experiment wo therefore need to find 
the area which gives the highest SNR given a set of survey 
and instrument parameters. 

This optimization procedure could be done for amy cos- 
mological parameter, or combination of parameters. How- 
ever for simplicity, and as a prerequisite to the measure- 
ment of the polarization power spectra, we will maximize 
the SNR for the amplitude. In principle, other parameters 
for the survey or telescope could be left free, such as the pix- 
els size or beam width. In practice we find that the smallest 
pixel/beam size is preferred, and so we set this to the limit 
of a given experiment. 



6.2 QUaD instrument parameters 

As a specific example of a ground-based experiment we use 
the QUaD experiment. This enables us to fix the instrument 
parameters needed to determine the pixel noise (equation 
13) and the allowed multipole range. These parameters are 
given in Table 2. A detailed description of QUaD is given in 
Church et al. (2003). 

The maximum multipole which can be covered is lim- 
ited by the beamsize as discussed in Section 3. If no other 
effect needs to be taken into consideration the minimum 
multipole, £min, would be determined from the survey area, 
imin = tt/G. Howcver, for a ground-based experiment the 
lower-£ cut-off is also limited by the stability of the atmo- 
sphere. This will limit the maximum scan which can be used 



Frequency (GHz) 


100 


150 


Number of bolometers 


21 


38 


Angular resolution (arcmin) 


6.3 


4.2 


NET per bolometer'^ {^lKs^^^) 


270 


300 



a - The definition of sensitvity for a polarization sensitive bolome- 
ter is discussed in Appendix A. 



and hence the largest angle on the sky over which a correla- 
tion can be made. For a perfect polarization experiment, this 
would not be an issue, as the unpolarized atmospheric fluctu- 
aions would not be detected in the polarized data. However, 
instrumental effects will cause a fraction of the unpolarized 
(common mode) signal to be to be present in the polarized 
signal. The atmosphere at the South Pole is exceptionally 
stable (Lay & Halverson 1998) and the QUaD instrument 
has been designed is such a way that these effects will be 
minimized, so we estimate that a minimum i of 25 can be 
reached. The minimum £ used in equation (28) will then 
be imin = max(7r/0, 25). Although it is possible for QUaD 
to make total power measurements, there is no mechanism 
for removing the atmospheric noise from the resulting data 
and so we assume that QUaD would not be able to produce 
temperature maps. 

We estimate the total observing time by assuming that 
QUaD will observe at the South Pole for two years during 
the austral winter (six months per year) for 22 hours each 
day and assuming that 20 per cent of this total time will be 
lost due to bad weather, instrument maintenance and cali- 
bration time. These estimates are based on the experiences 
of the DASI team at the South Pole site (Kovac et al. 2002). 
This gives a total time spent observing on the CMB of 3210 
hours per year. The maximum useable patch of sky is about 
1000 deg^, limited by available sky visible from the survey 
site, and major foreground contamination from the Galactic 
plane (see Fig. 3). We therefore restrict the analysis to areas 
below this maximum survey size. 

It is important to note that to measure both the Q and 
U Stokes parameters, each pixel must be meeisured with the 
detector in at least two different orientations with respect 
to the sky. For QUaD this will be achieved by rotating a 
half-wave plate so that both Q and U can be measured by 
each detector. This halves the total integration time avail- 
able for each Stokes parameter when making a polarized 
measurement. 



6.3 Results 

We have applied the above procedure to a model QUaD 
experiment. We consider three different cases: 

(i) a measurement of the EE spectrum, 

(ii) a measurement of the BB spectrum including the lens- 
ing component as part of the signal we wish to measure, 

(iii) a measurement of just the BB GW spectrum includ- 
ing the lensing signal as an extra source of noise. 

The results for the QUaD parameters are presented in 
Table 3, for TE, EE and BB spectra, for the case of fore- 
grounds, and without foregrounds. The latter is of interest 
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Table 3. SNRs for the optimal survey areas for each of the power 
spectra for a two-year integration time with QUaD 
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Figure 4. Variation of SNR with survey area for the EE (left) 
and total BB (right) power spectra with foregrounds (solid lines) 
and without foregrounds (dashed lines) for an observing time of 
1 year (light) and 2 years (dark) for r=0.01 



if the foregrounds are well enough understood to be sub- 
tracted from the signal or if a patch of sky with very low 
foreground variance can be found, as discussed in Section 4 
Fig. 4 shows how the SNR varies with area for EE and 
BB spectra. From Fig. 4 (left) it can be seen that a ground- 
based polarization experiment can make a good measure- 
ment of the EE spectrum, even if the foreground contami- 
nation is not well understood. The SNR is close to 100 for 
survey areas over 300 dcg^. Below this size the SNR falls 
rapidly to zero. For an experiment with the sensitivity and 
multipole coverage of QUaD, an E-mode survey is sample- 
variance-limited, so that a larger area is preferable (> 1000 
deg^) for statistical purposes. Given that an E-survey will 
have high signal-to-noisc per pixel, a high-resolution polar- 
ization map of the surveyed area is also possible (see Section 
7), allowing for removal of point sources as discussed in Sec- 
tion 4. 

Fig. 4 (right) shows the SNR for a two-year B-mode 
survey. The dotted Hue is for a SNR=3 which is the min- 
imum SNR that can be considered as an actual detection 
of the signal. Unlike the El-mode survey, the B-mode sur- 
vey is detector noise limited as the signal is much lower. 
With foregrounds the SNR sharply peaks at SNR ~ 5.5 for 
much smaller areas, around 50 deg^, where the lensing sig- 
nal dominates. As the survey area increases the noise per 
pixel increases and the overall SNR drops. If we can remove 
foreground contamination then the maximum SNR increases 
to a value of 9 and the optimal area is slightly reduced, as 
shown in Table 3. 

This different behaviour between the E and B-mode sur- 
veys with increasing area makes the simultaneous optimiza- 
tion of both measurements difficult. One compromise is to 
use the break in SNR of the E-mode survey at around 300 
deg^ . Although this is sub-optimal for both surveys, the drop 
to SNR = 90 for the E-survey is minimal and SNR = 4.5 for 
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a - including astrophysical foregrounds 
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the B-mode survey is still a strong detection. An alternative 
would be to split the survey in two, one large, one small, 
halving the integration time for each survey. As shown in 
Fig. 4 for a single year of integration it is still possible to 
detect the B-mode signal if we concentrate on a small area 
of sky (~ 20 deg^). For a single year the EE SNR also does 
not drop significantly if an area larger than about 500 deg^ 
is chosen. 

From Table 3 it is evident that QUaD cannot detect the 
GW B-mode component unless the tensor to scalar ratio is 
larger than the values considered so far. We have therefore 
extended the calculation to higher values of r up to the cur- 
rent upper limit. Fig. 5 shows how the optimal area for a 
measurement of the GW signal with QUaD varies with r. 
The optimal area changes significantly as r increases. For 
the foreground model assumed here it is only possible to de- 
tect the GW signal for r greater than 0.35. However, for the 
large areas which are best for detecting this high GW signal, 
the SNR for the total B-mode signal drops significantly. It 
is therefore not possible to pursue both science goals simul- 
taneously. However, if the foreground comtamination can be 
completely removed, the lowest detectable value of r drops 
to 0.14. The optimal area also decreases as the detector noise 
becomes the dominant factor. In the no foreground case it 
would be possible to detect the GW signal using the 300 
deg^ survey discussed above. If the GL signal can be re- 
moved, the GW signal becomes slightly easier to detect, but 
only if the foregrounds can be subtracted as the combined 
dust and synchroton contamination (Fig. 2) is laxger than 
the GL signal over most of the multipole range which can 
be covered from the ground. 

The effects of the mixing of E and B modes due to par- 
tial sky coverage will not significantly influence the results 
found here. Bunn (2002) finds that the mixing will only have 
a large effect for the B-mode signal on the scale of the sur- 
vey size. If we use a 300 deg'^ patch the GL B-mode signal 
will therefore not be affected. For a detection of the GW sig- 
nal this effect will become more important. However, Lewis, 
Challinor & Turok (2002) discuss this problem and calculate 
the minimum detectable r as a function of survey size. They 
find that for the large surveys (greater than 50°) the mini- 
mum value is not changed if the mixing effects are included. 
For the areas discussed here the GW results will therefore 
not be influenced by E-B mixing if an optimal method is 
used to separate the E and B modes. 

To estimate the SNR for the TE spectrum we assume 
that a QUaD map could be combined with the portion of 
the expected four-year WMAP data covering the same area 
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Table 4. Planck and WMAP instrument parameters. For WMAP 
we use only the highest two frequency channels as the other chan- 
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Figure 5. Variation of the survey optimal area (upper plot) and 
achievable SNR (lower plot) for a measurement of the GW sig- 
nal with QUaD as a function of r with (solid line) and without 
(dashed line) foregrounds. 



of sky. The results axe shown in Table 3. As with the EE 
spectrum, the measurement is sample-variance-limited and 
the largest possible area of 1000 dcg'^ is best. The SNR also 
drops sharply if the survey area becomes too small (< 100 
deg^). However, the QUaD TE measurement is limited by 
the resolution and sensitivity of the WMAP map and suffers 
more heavily from sample variance than the smaller EE sig- 
nal. The SNR with which this spectrum could bo measured 
by QUaD is therefore smaller than the EE SNR. The TE 
spectrum has also already been measured in this multipole 
range by WMAP. It is therefore more useful to optimize a 
ground-based survey for a measurement of the EE and BB 
spectra. 

We have also investigated the effect of increasing the 
minimum I value used in the calculation. For the TE, EE 
and total BB spectra, an increase in the minimum I from 25 
to 100 has a negligible effect, as most of the power in these 
spectra is from the higher multipoles. However, as would bo 
expected, increasing the minimum i does affect the GW B- 
mode detection. If the minimum I is increased to 100 the 
GW is no longer detectable below the current upper limit of 
0.36. 



7 DEEP MAPS OF THE CMB POLARIZATION 

For a ground-based experiment it is not possible to make ob- 
servations of the whole sky due to the limited sky coverage 
available from the ground. Although this is a disadvantage 
in terms of multipole coverage at low I, by making a deep in- 
tegration of a small region of sky it is possible to make maps 
with a very high signal-to-noise ratio. This allows more pre- 
cise measurements to be made on small angular scales. It 
will also improve the ability of the experiment to remove 
low-lying systematic effects which would not be detectable 
in observations of lower signal-to-noise. We illustrate the dif- 
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ference between QUaD and the Pland^^ satelUite mission 
in Fig. 6 using simple simulations of the Q Stokes parame- 
ter with noise appropriate to each experiment. While Planck 
will cover much more sky than QUaD, the QUaD observa- 
tions would be at higher signal-to-noise than the average of 

those made by Planck^^ 1 1 . This will allow QUaD to limit sys- 
tematic effects in the experiment to an unprecented level. 
This makes the QUaD approach highly complementary to 
that of Planck, which would have lower average sensitiv- 
ity, but good statistics over the entire sky. The deep maps 
will also provide new information on the techology used by 
QUaD and on polarized foregrounds, which will be crucial 
for the design of future CMB experiments. 

Using a smaller region of sky is also an advantage in 
terms of foregrounds, as it is possible to target the most 
useful patches of sky, without spending valuable integration 
time on regions which will ultimately be left unused in cos- 
mological analyses. 

Finally, it is possible tailor the size of the region ob- 
served to optimize for a particular science goal. As was de- 
scribed in Section 6, this is especially important for searches 
for the faint B-mode signal. 



8 POWER SPECTRA ESTIMATION 

By statistically averaging over the polarization signal, the 
polarization power spectra may be estimated. Fig. 7 com- 
pares the expected band-averaged power spectra results and 
multipole coverage of a 300 deg^ , two-year ground-based ex- 
periment, QUaD, an all-sky four-year satellite, WMAP, and 
the Planck satellite mission. These predictions are based on 
equation (11), with the parameter the power in a pass-band 
of width A£. We include all of the power covariances and 
effects of foreground emission outlined in Section 4. The in- 
strument parameters used in each experiment are given in 
Tables 2 and 4. 

From this analysis we find that QUaD can make a high- 
significance measurement of the EE-power spectrum, as sug- 
gested by the high SNR (~ 90) found during optimization 
in Section 6, over a multipole range from / = 25 to = 2500. 
The polarization acoustic oscillations are well sampled, with 
a resolution of A£ w 20. The lower modes are not sampled 



11 http: / / www. astro, esa. int / S A-general/Proj ccts /Planck / 
III! It is noted that the average noise over the entire sky was used 
for Planck. While Planck will cover some regions, namely the 
Ecliptic poles, more deeply, these regions are in general not the 
best in terms of foregrounds, and the mean noise away from these 
regions will be correspondingly worse 
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Figure 6. Simulated polarization maps from QUaD (left) and Planck (right), showing a 300 deg^ field of Q-mode anisotropies at 150 
GHz. All of the structure in the high signal-to-noise QUaD map is real signal, while the Planck map (shown with the same pixelization) 
has much lower signal-to-noise per pixel. Foregrounds and other systematic effects are not included in either map. 
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Figure 7. Predicted measurements of the polarization power spectra achievable with the current generation of satellite (four-year WMAP 
and two-year Planck) and ground-based (two-year QUaD) experiments for r = 0.1. The error bars show detections above the one-sigma 
level and the free symbols show upper limits. The errors include a contribution from astrophysical foregrounds. For clarity we do not 
show the Planck measurements of the EE power spectrum. These will be of a similar sensitivity to the QUaD measurements, but will 
cover different I ranges, as indicated by the points on the B-mode spectrum. 



due to the limited survey area. In particular the re-ionization 
peak at £ = 7 is only detected by a satellite mission. 

In addition there is a good detection of the BB-power 
spectrum from £ = 25 to ^ = 1000. Power is binned loga- 
rithmically to increase the signal-to-noise per bin. The most 
significant bin is at f = 1000, at the peak of the gravita- 
tional lensing (GL) contribution to the BB-spectrum. If the 
foreground contamination can be significantly reduced, a di- 
rect detection of the GW contribution to the B-mode power 



spectrum could be made at around I = 100. Again the low-^ 
modes are not accessible to a ground-based survey, but can 
be complementarily detected by an all-sky satellite mission. 

With both temperature and polarization data available 
the TE-cross power spectra may also be estimated so that 
a cross-check can be made with other measurements of this 
signal. 

With such high-resolution polarization information 
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available it is interesting to see what effect a ground-based 
survey will have on cosmological parameters. 



9 PARAMETER ESTIMATION 

In this Section we investigate the contribution which can 
be made by ground-based polarization experiments to the 
measurement of the cosmological parameters. Previous work 
on CMB parameter estimation (Efstathiou & Bond 1999; 
Zaldaxriaga, Spergel & Seljak 1997; Bond, Efstathiou & 
Tegmark 1997) has shown that the polarization data which 
can by obtained by the forthcoming WMAP and Planck 
satellite missions will allow a more accurate determination 
of many of the key cosmological parameters. For a satellite 
experiment, this is mainly because the degeneracy between 
T and A can be broken by measuring the re-ionization bump 
in the polarization power spectra. These re-ionization bumps 
also create a high GW B-mode signal at low £, so a full-sky 
measurement will also tighten constraints on the tensor-to- 
scalar ratio. 

As we have discussed, a ground-based polarization ex- 
periment can concentrate on smaller areas of sky at higher 
resolution and so can make a good measurement of the 
acoustic peaks out to high £ in the EE-power spectrum. The 
information from a ground-based experiment will therefore 
complement the full-sky satellite data. It is also possible 
to choose an observing strategy with targets the GW sig- 
nal peak at intermediate scales {i = 100). Ground-based 
constraints on the B-mode GW signal will therefore also 
complement those obtainable from the current generation 
of satellite experiments. 

Finally, it is important to note that a CMB polarization 
experiment is not just adding more data. A similar exper- 
iment measuring only the temperature spectrum, over the 
same multipole range, and with the same detector sensitiv- 
ity, would add very little new information as far as cosmo- 
logical parameters are concerned, although a high-£ tem- 
perature surveys may well start to probe higher-order CMB 
effects. Hence polarization adds unique information from the 
CMB. 

We investigate the potential increase in the precision of 
the measurement of cosmological parameters which can be 
achieved with a ground-based experiment by comparing the 
expected four- year results from WMAP alone to those which 
could be achieved by combining QUaD and WMAP data. To 
compare the two cases we calculate the inverse Fisher ma- 
trix using equation (9) to find the variances and covariances 
between each of the parameters. For QUaD we use the in- 
strument model dicussed in Section 6.2. The experimental 
parameters used for WMAP are given in Table 4. 

To calculate the derivatives in parameter space required 
in equation (9) we use second-order differencing between 
CMBFAST models for accuracy, with the corresponding pa- 
rameter changed up and down by 1 per cent. The deriva- 
tive of the power spectrum encapsulates the response of the 
spectrum to a change in a particulax parameter and hence 
quantifies its information content. However, if the shape of 
the derivative for any two parameters is too similar then the 
two parameters will be degenerate and cannot both be con- 
strained. The derivatives used in the calculation are shown 
in Fig. 8 and 9. For most parameters, the shape of the deriva- 



tives reflects the acoustic peaks in the power spectrum, in- 
dicating that both information about parameters and their 
differences are contained in the peaks. For instance with 
temperature only Qmh'^ and Os/i^ are quite anti-correlated, 
but their derivatives oscillate out of phase for TE and EE- 
spectra, breaking this degeneracy. Much of the difference 
between h and ris occurs in the low multipoles, but there 
is a large difference at high-£ in the BB-spectra due to the 
effects of gravitational lensing on these modes. 

Fig. 9 clearly shows the anti-correlation that arises be- 
tween T and A when only temperature information is avail- 
able. This degeneracy is seen to be broken on large scales 
by the differences in the responses of the polarization power 
spectra. However, going to high £ in the polarization spec- 
tra, these parameters become strongly degenerate again. 
Hence we can expect that a ground-based polarization sur- 
vey, which will have difficulty reaching the lower multipole 
range, will not contribute much to lifting the A — t degener- 
acy. Conversely, with only temperature information, h and 
ris are strongly degenerate, with much of the difference in 
response coming at very low modes, or modes beyond a few 
hundred. However, adding polarization information, espe- 
cially TE at around £ of 100, and lerised BB modes at high 
£, breaks this degeneracy due to their different responses 

Finally with only a temperature spectrum, the response 
to r is limited to the first hundred multipoles. The TE and 
EE derivatives show that there is useful information about 
r on intermediate scales in the polarization spectra up to 
around £ = 500, but for these multipoles the scalar EE and 
TE power spectra are very high and so it will be difficult to 
extract this information from the signal. However, the BB 
derivative also shows structure at higher multipoles and will 
provide information on r if the tensor signal is higher than 
the scalar lensing signal at the scales of interest, or if the 
lensing signal can be removed. 

Having considered the responses of the power spectra to 
our parameter set, we now turn to estimating parameter un- 
certainties from satellite and ground-based surveys. To test 
the validity of this procedure we have calculated the accu- 
racy achievable with the one-year WMAP data, and found 
that our results are on good agreement with the one-year 
WMAP quoted parameter errors (Spergel et al. 2003). 

Fig. 10 shows the relative error ellipses (defined by the 
AlnL = —1/2 contour) expected from a 4-year WMAP 
experiment (darker ellipses) and from a combined 2-year 
ground-based QUaD and 4-year WMAP experiment (lighter 
ellipses) for our fiducial 7-parameter model, marginalizing 
over the other parameters. The projection of this contour 
gives the marginalized one-parameter, 1 — a error for each 
parameter. For a two-parameter 68 per cent confidence re- 
gion, the ellispes should be scaled by a factor 1.5. We assume 
that the TE-cross spectrum can be estimated in the overlap 
region. Here we see that a significant improvement of around 
a factor 2 is made on most of the parameter set by adding in 
a ground-based polarization survey, despite the significant 
difference in survey size. For most parameters, this comes 
from the high-multipole information in the EE-spectra, but 
there is also important information in the BB-spectra, in 
particular for r, h and Ug. 

The 1-cr marginalized parameter uncertainties for 
WMAP and QUaD -|- WMAP are shown in Table 5. By 
including QUaD the precision with which the parameters 
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Figure 8. The derivative of the CMB polarization power spectra with the parameters Hg/i^, Clrnh^, h, and ris- The models are generated 
by CMBfast, with fiducial parameters given in the text. 
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Figure 9. The derivative of the CMB polarization power spectra with the parameters r, A and r. The models are generated by CMBfast, 
with fiducial parameters given in the text. Note that the plot is now logarithmic in order that the low i degeneracy breaking between r 
and A can be observed in the polarization power spectra. 



can be mesisured is improved by around a factor of two in 
most cases. This increase in accuracy arises from the ex- 
tra information in the EE-spectra from modes £ > 100, and 
from the strong BB-spectral dependence on small scales for 
ris and h. Again, for a temperature survey alone Fig. 8 and 
9 indicate there is no useful information at high-multipoles. 



It is interesting to look at how the information from 
the B-mode spectrum influences the parameter estimation. 
To examine this the same calcuation was made, but with 
the B-modc information removed from the Fisher matrix. 
For WMAP, this did not change the parameter estimates 
significantly, except for a slight increase in the error on r 
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Figure 10. Marginalized Fisher matrix relative parameter error constraints (AlnL = —1/2) anticipated for four-year WMAP results 
only (dark) and four-year WMAP combined with QUaD (light) for r = 0.01 with foregrounds. The projections of the ellipse onto the 
two axes give the standard errors on each parameter. For a two-parameter 68 per cent confidence region, the ellipses should be scaled by 
a factor 1.5. 



(~ 10 per cent). For WMAP most of the information on r 
must therefore come from the TT, TE and EE spectra, and 
not from the weak upper limit on the B-mode spectrum. 
For QUaD we find a shght increase in the errors on h and 
Us (~ 20 per cent) due to the loss of the information con- 
tained in the B-mode lensing signal. However, the error on r 
more than doubles if B-modes are not included. The B-mode 
information from QUaD must therefore make a significant 
contribution to the r constraint, even though QUaD cannot 
make a strong detection of the GW B-mode signal. 

Fig. 11 shows the predicted improvement on a joint 
measurement of r and Us from a two-year QUaD experi- 
ment and four-year WMAP survey. With a detection of r, 
Us and the amplitude A, the shape of the inflaton potential 
can be inferred (Hoffman & Turner 2001). 

The poorest parameter improvement is for t and A, 
which only improve by a factor of about 1.3. As discussed 



above this is because the main differences appear on scales 
oi i < 100, which are difficult to reach from the ground, but 
accessible to satellite surveys. 



We find that for most parameters the errors do not de- 
crease significantly if the foreground contamination is com- 
pletely removed. However, this is not the case for r, where 
the error decreases by a factor of three if the foreground 
contamination can be removed, leading to a factor of six 
improvement over the WMAP-only constraints. This can be 
clearly seen from the inner contours in Fig. 11. This is be- 
cause the foreground removal allows a much better measure- 
ment of the B-mode GW signal to be made. 
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Table 5. Fisher matrix estimates of parameter errors 



Parameter 


Value 


WMAP 


WMAP +QUaD 




0.0224 


0.0009 


0.0004 




0.135 


0.007 


0.004 


h 


0.71 


0.040 


0.021 


T 


0.17 


0.023 


0.020 


Us 


0.93 


0.029 


0.014 


r 


0.01 (0.1) 


0.206 (0.203) 


0.082 (0.090) 


A 


0.83 


0.036 


0.031 






CO 







Figure 11. Predicted improvment on a joint measurement of r 

and Us for two-year WMAP data, (outer contour), two-year QUaD 
+ WMAP including foregrounds (middle contour) and two-year 
QUaD + lyA'MP without foregrounds (centre contour) for r = 0.1 
(left) and r = 0.01 (right). The AlnL = —1/2 contour is shown. 
For 68 per cent confidence limits, scale by a factor 1.5 



10 SUMMARY 

In this paper we have investigated the science goals achiev- 
able with the forthcoming generation of ground-based CMB 
polaxization experiments. We have set out a Fisher infor- 
mation matrix formalism that takes into account the com- 
bination of different temperature and polarization surveys, 
and includes foreground contamination. We have argued 
that ground-based polarization experiments can reach the 
high sensitivities required by making a deep integration on 
a small patch of sky. By preferentially selecting regions of 
sky with low foreground variance it will also be possible for 
a ground-based experiment to reduce further the foreground 
contamination. 

Taking the proposed QUaD South Pole experiment as 
our model survey, we have optimized the survey area and 
shown that a 300 deg^ survey is a good compromise be- 
tween a sample-limited E-mode survey and a detector-noise 
limited B-mode survey. Below 300 deg^ the SNR for the E- 
mode survey drops rapidly, while above this a detection of 
the (gravitational lensing component of the) BB-power spec- 
trum becomes unfeasable. With such high signal-to-noise 
per pixel in the E-mode survey, deep imaging maps of the 
CMB polaxization field can be made. Statistically averaging 
the data allows a high-significance measurement of the EE- 
power spectrum over a range of multipoles from £ = 25 to 
i = 2500, with good sampling of the acoustic oscillations. 
The gravitational lensed component of the BB-power spec- 
trum can also be detected with good signal-to-noise. If it is 



possible to reduce the foreground contamination the gravi- 
tational wave could also be detected for r > 0.14. 

Combining a two-year QUaD experiment with a four- 
year WMAP all-sky survey allows a better measurement of 
cosmological parameters to be made compared to that pos- 
sible from WMAP data alone. Most parameters can be im- 
proved by a factor two. If the foreground contamination can 
be reduced the tensor-to-scalar ratio will be dramatically 
improved by up to a factor of six. With such improvements, 
strong constraints can be placed on the potential of the in- 
flaton field. Only the degeneracy between the amplitude of 
fluctuations, A, and the optical depth to re-ionization, t, are 
not significantly improved, as this requires large scales only 
accessible to a satellite. 

In conclusion we find that if the necessary sensitivity 
and control of systematics can be achieved, a ground-based 
CMB polarization experiment such as QUaD can make a 
major contribution to the study of CMB polarization power 
spectra and cosmological parameters. 
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APPENDIX A: SENSITIVITY DEFINITIONS 
FOR CMB POLARIZATION EXPERIMENTS 

There are a number of definitions for the sensitivity of a 
CMB polarization experiment and this is often a cause of 
confusion when comparing different sensitivity parameters. 
For a total- power CMB experiment the sensitivity is usually 
defined in terms of the noise equivalent temperature (NET) . 
This is the signal needed from the source to give a signal- 
to-noise ratio of unity in a one-second integration time.* * * 
To measure polarization an equivalent definition is required 
in terms of the Q and U Stokes parameters. For a linearly 
polarized source of total intensity, 7, of which a fraction p is 
polarized at an angle x to the reference direction, the Stokes 
parameters can be defined as: 



Q =p/cos(2x), 



* * * For CMB work the sensitivity is usually quoted as an NET, 
in units of Ks^^^, instead of as a noise equivalent power (NEP), 
which is normally used in sub-millimetre astronomy. This makes 
it easier to combine experimental work with theory, as the power 
spectra (C^) are defined in terms of temperature units. The 
NEP is normally quoted per unit bandwidth and so has units of 
WHz~^^'^, which is equivalent to noise produced in a half second 
integration time. To change NET in ifsV2 to NEP in WRz'^/'^ 
the conversion is: 



NET{Ks^/'^) : 



NEP{WHz-'^/^) 

V2 dB^/dT ' 



(Al) 



where dB^/dT is the derivative of the source (the CMB) with 
respect to temperature. The factor of \/2 converts from Hz to 
seconds. 
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[/ =pJsin(2x). 



(A2) 



NEQ = V2vdNET, = 



If we orientate the axis of the reference system so that it 
is aligned with the polarization angle of the source (x = 0) 
then we have U = and Q = pi so that Q gives the total 
polarized intensity. We can then define the polarization sen- 
sitivity, NEQ, as the polarized signal from the source needed 
to give a signal-to-noise ratio of unity in a one-second inte- 
gration time for a source with a polarization angle aligned 
with the reference direction of the measurement. 

For QUaD, the polarized measurements will be made 
with pairs of polarization-sensitive bolometers (PSBs). The 
two bolometers in a PSB pair arc sensitive to orthogonal 
polarization states of the incoming radiation. The intensity 
measured by the co-polar (a;) and cross-polar (y) device is 
given by (Jones et al. 2003): 



(A8) 



Q), 



(A3) 



in a reference system aligned with the polarization angle of 
the source. The total intensity is found by adding the two 
bolometer outputs and the Q Stokes parameter is found by 
differencing the outputs. 

For a bolometer, the noise equivalent power due to pho- 
ton noise (NEP) is given by (Lamaxre 1986): 

(A4) 

where is m is the number of polarization states detected (m 
is either 1 or 2). For a single PSB, m = 1, as only a single 
polarization state is detected. P is the power in a band of 
width Aj/: 



9 P 

NEP^ = 2hvP + 



Pv = ridritsAflBr^Au, 



(A5) 



where Afl is the throughput of the system, e is the emis- 
sivity of the source and is the intensity of the radiation 
that would be emitted from a perfect black body. The total 
efficiency of the system is = ridrjt , where ija is the detector 
efficiency and rji is the instrument efficiency. We assume that 
a lossless PSB will absorb half of the incident unpolarized 
radiation, giving r)d = 1/2. The NET due to photon noise 
in each PSB from the unpolarized background radiation is 
therefore: 

^^^^ = ^ ,f,' . A6 

Tjdrit OBv/aT 

where the factor of ridfit is needed to convert from the noise 
at the detector to the signal required at the source, as only 
a fraction of the radation from the source will be absorbed 
by the PSB. dBu/dT is the derivative of the source intensity 
with respect to temperature and converts from an NEP to 
an NET. Equation (A6) gives the NET for a measurement 
of the temperature of the CMB with a single PSB. 

In order to measure the polarization we require a pair 
of PSBs. The temperature sensitivity of a PSB pair can be 
obtained by averaging the two outputs so that: 

NETj,air = NETs/V2. (A7) 

NETpair is exactly the NET that would be obtained if a sin- 
gle normal (not polarization sensitive) bolometer had been 
used. For a measurement of Q, the two outputs are differ- 
enced so that: 



An important point to note is the factor of rid in this expres- 
sion. This is because we are now measuring the signal from 
a polarized source, so the factor of rja which was needed in 
equation A6 to find the noise for a measurement of the total 
power is no-longer required. All of the polarized radiation is 
absorbed by a single PSB when it is correctly aligned with 
the polarization angle of the source. 

When defining the sensitivity of a PSB it is therefore 
important to state whether a sensitivity is an NET for a 
single detector, an NET for a pair of detectors, or an NEQ 
for a pair of detectors. The expression for the pixel noise 
given in Section 3 will depend on the sensitivity definition 
used: 

2 NET^Q'^ ,^ , 

tobsJypsBf'ip,^ 

If the NET is for a single PSB (as in Table 2 ), Npsb is the 
total number of PSBs. If the sensitivity is given as an NEQ 
for a PSB pair, Npsb is the number of pairs. 
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